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Abstract 

We determine tlie hard-loop resummed propagator in an anisotropic QCD plasma in general 
covariant gauges and define a potential between heavy quarks from the Fourier transform of its 
static limit. We find that there is stronger attraction on distance scales on the order of the inverse 
Debye mass for quark pairs aligned along the direction of anisotropy than for transverse alignment. 
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I. INTRODUCTION 



Information on quarkonium spectral functions at high temperature has started to emerge from 
lattice-QCD simulations; we refer to ref. [l| for recent work and for links to earlier studies. This has 
motivated a number of attempts to understand the lattice measurements within non-relativistic 
potential models including finite temperature effects such as screening A detailed discussion 
of the properties of the heavy-quark potential in the deconfined phase of QCD is given in ref. Q], 
which also provides a comprehensive list of earlier work. Also, Laine et al. have recently derived a 
Schrodinger equation for the finite-temperature Wilson loop to leading order within "hard-thermal 
loop" (HTL) resummed perturbation theory by analytic continuation to real time 0]. Aside from 
the well-known screened Debye potential, their result includes an imaginary part due to Landau 
damping of low-frequency modes of the gauge field, corresponding to a finite life-time of quarkonium 
states. 

The present paper is a first attempt to consider the effects due to a local anisotropy of the 
plasma in momentum space on the heavy-quark potential. Such deviations from perfect isotropy are 
expected for a real plasma created in high-energy heavy-ion collisions, which undergoes expansion. 
The HTL propagator of an anisotropic plasma has been calculated in time-axial gauge in ref. 0]. We 
derive the result for general covariant gauges, which allows us to define a non-relativistic potential 
via the Fourier transform of the propagator in the static limit. 



II. HARD-THERMAL-LOOP SELF-ENERGY IN AN ANISOTROPIC PLASMA 



The retarded gauge-field self-energy in the hard-loop approximation is given by 



p - V + ie 



(1) 



Here, = (l,k/|k|) is a light-like vector describing the propagation of a plasma particle in space- 
time. The self-energy is symmetric, Yl^^^{p) = n^'^(p), and transverse, p^Yi^^{p) = 0. 

In a suitable tensor basis the components of H^'^ can be determined explicitly. For anisotropic 
systems there are more independent projectors than for the standard equilibrium case [sl. Here, 
we extend the tensor basis used in to a four-tensor basis appropriate for use in general covariant 
gauges. Specifically, 
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Here, is the heat-bath vector, which in the local rest frame is given by = (1,0,0,0), and 

= - p^" (6) 
p^ 

is the part that is orthogonal to p^ . 
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The direction of anisotropy in momentum space is determined by the vector 

n'^ = (0,n), (7) 

where n is a three-dimensional unit vector. As before, is the part of orthogonal to p^. 
The self-energy can now be written as 

U'"' = aAi"^ + 133^" + 7C^^ + dDf"" . (8) 

In order to determine the four structure functions explicitly we need to specify the phase-space 
distribution function. We employ the following ansatz: 



/(p) = /iso(^VP' + ^(P-")'j • (9) 
Thus, /(p) is obtained from an isotropic distribution /iso(|p|) by removing particles with a large 
momentum component along n. The function /iso(|p|) should decrease monotonically with |p|, so 
that the square of the Debye mass defined in eq. below is guaranteed to be positive; however, 
in the real-time approach employed here, the distribution /iso need not necessarily be thermal. 

The parameter ^ determines the degree of anisotropy, ^ = (1/2) (pj^)/ (pi) — 1, where = p • n 
and p± = p — n(p • n) denote the particle momentum along and perpendicular to the direction n 
of anisotropy, respectively. // /iso is a thermal ideal-gas distribution and ^ is small then ^ is also 
related to the shear viscosity of the plasma; for example, for one-dimensional Bjorken expansion 0] 

« = |^5. (10) 

Its 

where T is the temperature, r is proper time, and rj/s is the ratio of shear viscosity to entropy 
density. In an expanding system, non- vanishing viscosity implies finite momentum relaxation rate 
and therefore an anisotropy of the particle momenta. 

Since the self-energy tensor is symmetric and transverse, not all of its components are indepen- 
dent. We can therefore restrict our considerations to the spatial part of 11'^'^, 

{p, = ml -r^ jT—w 4^ + T- ' ^1 

and employ the following contractions: 

^^'n^nV = (P^ - {n ■ pf)5 , 
P^ 

TYff^' = 2Q + /3 + 7 . (12) 
The Debye mass mo appearing in eq. (jlip is given by 

2 9^ r , 2dhM 



^n = -^,l dpp (13) 

where p = |p|. We do not list the rather cumbersome explicit expressions for the four structure 
functions a, /3, 7, and 5 here since they have already been determined in ref. [H]- 

In principle, the tensor basis dSJlS]) could be chosen differently, such that the individual tensors 
have a simpler structure. For example, one could choose 

Tn • 71 

C^"" = h^n'' - ^^{m}'n'' ^m'^n^) (14) 

j^i^v = _ ^M^^ (15) 

However, in the basis ([2][5|) the spatial components of 11'^'' are identical to those from ref. [5] and 
so we can avoid the rather tedious re-evaluation of the four structure functions. 
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III. PROPAGATOR IN COVARIANT GAUGE IN AN ANISOTROPIC PLASMA 



From the above result for the gluon self-energy one can obtain the propagator iA[^^ . It is 
diagonal in color and so color indices will be suppressed. In covariant gauge, its inverse is given by 



A 



(16) 



where co = p ■ m and A is the gauge parameter. Upon inversion, the propagator is written as 

A'^^(p, ^) = a'Ai"" + 13' Bi"' + -f'Ci"' + + rjpi'p'' . (17) 

Using {A-^Y^A^ " = gi"" it follows that the coefficient of gf"" in (A-^)'''"A^ " should equal 1 while 

the coefficients of the other tensor structures, for example of n^n^ , n^p^ and p^p'^, should vanish. 
Hence, we can determine the coefficients in the propagator from the following equations 

1 
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Hence, we find that in covariant gauge the propagator in an anisotropic plasma is given by 
1 
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For ^ = 0, we recover the isotropic propagator in covariant gauge 
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IV. HEAVY QUARK POTENTIAL IN AN ANISOTROPIC PLASMA 

We determine the real part of the heavy-quark potential in the nonrelativistic limit, at leading 
order, from the Fourier transform of the static gluon propagator, 

V{r, = -g'CF J e^P-^AOO(u; = 0, p, (26) 
(2^0^ ''''' ' (P' + ml + m2)(p2 + In}) - mj ' ^^^^ 
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The masses are given by 
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and 



P'=pi+P'- (32) 

The above expressions apply when n = (0, 0, 1) points along the 2:-axis; in the general case, Pz and 
p_L get replaced by p • n and p — n(p • n), respectively. 

We first check some limiting cases. When ^ = then mp = ruD while all other mass scales in 
the static propagator vanish. Hence, we recover the isotropic Debye potential 

nr,^ = 0) = V,M = -9C,J ^:^^ = -^e , (33) 

where r = rmn- 

Consider, on the other hand, the limit r — > for arbitrary ^. The phase factor in ()27p is 
essentially constant up to momenta of order |p| ~ 1/r and since the masses are bounded as 
IpI — > CO they can be neglected. The potential then coincides with the vacuum Coulomb potential 

V{r ^ 0,0 = Kac(r) = -g^C, J = . (34) 

The same potential emerges for extreme anisotropy since all mj — > as ^ — > oo: 

y(r,e = oo) = -^. (35) 

This is due to the fact that at ^ = oo the phase space density /(p) from eq. ([9]) has support only in 
a two-dimensional plane orthogonal to the direction n of anisotropy. As a consequence, the density 
of the medium vanishes in this limit. 

For an anisotropic distribution, the potential depends on the angle between r and n. This can 
be seen analytically for small but non-zero ^. To linear order in ^ the potential can be expressed 
as 

j^J^I/^'p ..„.J-(P-n)Vp^ 

For r parallel to the direction n of anisotropy, 



V{r, e « 1) = VUr) - g'^Cp iml / — ^ e^P- ^ , l\ ! " ■ (36) 
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(37) 



where r = rrriD, as before. This expression does not apply for f much larger than 1, which is 
a shortcoming of the direct Taylor expansion of V{r,^^) in powers of ^. However, for f ~ 1 the 
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FIG. 1: Heavy-quark potential at leading order as a function of distance (f = rm/j) for r parallel 
to the direction n of anisotropy. The anisotropy parameter of the plasma is denoted by ^. 
Left: the potential divided by the Debye mass and by the coupling, V = V/ (g'^CF'mo)- Right: 
potential relative to that in vacuum. 

coefficient of is positive, (• • •) = 0.27 for f = 1, and thus a slightly deeper potential than in an 
isotropic plasma emerges at distance scales r ~ l/mn- 
When r is perpendicular to n. 



The same limitations for f apply as in eq. (I37p . Here, too, the coefficient of the anisotropy parameter 
is positive, (•••)= 0.115 for f = 1, but smaller than for r || n. Hence, a quark-antiquark pair aligned 
along the direction of momentum anisotropy and separated by a distance r ~ l/m^i is expected 
to attract more strongly than a pair aligned in the transverse plane. 

For general ^ and f, the integral in (|27p has to be performed numerically. The poles of the 
function are integrable^ . In Fig. [1] we show the potential in the region f ~ 1 for various degrees 
of plasma anisotropy. One observes that in general screening is reduced, i.e. that the potential at 
1^ > is deeper and closer to the vacuum potential than for an isotropic medium. This is partly 
caused by the lower density of the anisotropic plasma. However, the effect is not uniform in the 
polar angle, as shown in Fig. [2j the angular dependence disappears more rapidly at small f, while 
at large r there is stronger binding for r parallel to the direction of anisotropy. Overall, one may 
therefore expect that quarkonium states whose wave- functions are sensitive to the regime f ~ 1 are 
bound more strongly in an anisotropic medium, in particular if the quark-antiquark pair is aligned 
along n. 

V. DISCUSSION AND OUTLOOK 

We have determined the HTL gluon propagator in an anisotropic (viscous) plasma in covariant 
gauge. Its Fourier transform at vanishing frequency defines a non-relativistic potential for static 



They are simple first-order poles which can be evaluated using a principal part prescription. 



V{r±n,C<^l) = Viso{r) l + C 




(38) 
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FIG. 2: Comparison of V{r \\ n,^^) and V{r _L n, ^). 



sources. We find that, generically, screening is weaker than in isotropic media and so the potential 
is closer to that in vacuum, in particular if the QQ pair is aligned along the direction of anisotropy. 

Our results are applicable when the momentum of the exchanged gluon is on the order of the 
Debye mass mo or higher, i.e. for distances on the order of \d = I/'tt-d or less. For realistic 
values of the coupling, ~ 0.3, Xd is approximately equal to the scale rmcd(^) ^ 0.5 (Tc/T) fm 
introduced in [1,0], where medium-induced effects appear. 

Following the discussion in ref. at short distances, r < r^e^^iT), the potential is given by 



where o" ~ 1 GeV/fm is the SU(3) string tension; color factors have been absorbed into the 
couplings. Since r^ed{T) ~ 1/T, it follows that at sufficiently high temperature rined(^) is smaller 
than y/a/a and so the perturbative Coulomb contribution dominates over the linear confining 
potential at the length scale Ad. Roughly, this holds for T ^ 2Tc. In this case, our result is 
directly relevant for quarkonium states with wavefunctions which are sensitive to the length scale 

On the other hand, for lower T the scale rmed(r) where medium- induced effects appear may 
grow larger than ~ y/a/a. In this regime, quarkonium states are either unaffected by the medium; 
namely, if the quark mass is very large and the typical momentum component in the wave function 
is l/rjnedCT")- Conversely, states with a root-mean square radius ^ r^cd{T) do experience 
medium modifications. For such states, however, it is insufficient to consider only the (screened) 
Coulomb-part of the potential which arises from one-gluon exchange. Rather, one should then 
sum the medium-dependent contributions due to one-gluon exchange and due to the string 
We postpone detailed numerical solutions of the Schrodinger equation in our anisotropic potential 
to the future. It will also be interesting to understand how the width of quarkonium states 0] 
which arises in HTL resummed perturbation theory due to Landau damping of modes with low 
frequency is affected by an anisotropy of the medium. 
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